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Abstract 

We study how the components of hadronic wave functions in light-cone quantization depend on 
the ultraviolet cut-off by relating them in a systematic way to the matrix elements of a class of 
quark-gluon operators between the QCD vacuum and the hadrons. From this, we derive an infinite 
set of scale-evolution equations for the individual contributions to parton distributions from the 
Fock expansion. When summed over all the contributions, we recover the well-known DGLAP 
equation. 
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In light-cone quantization and light-cone gauge, the hadronic states in QCD are expressed 
as an expansion of various quark and gluon Fock components [[I]. This expansion depends, 
among others, on the momentum cut-off used to truncate the theory, which is often inter- 
preted as the physical resolution scale. Although the physical observables, such as masses, 
angular momenta, form factors and cross sections, ought be independent of the cut-off, 
many interesting hadronic matrix elements do. A well-known example is the matrix ele- 
ments of twist-two operators which define the moments of Feynman's parton distributions 
0. For these quantities, one should be able to trace the scale dependence back to that of 
the hadronic wave functions. 

In this paper, we are interested in how the light-cone Fock components depend on the mo- 
mentum cut-off. Finding the solution directly from diagonalizing the light-cone hamiltonian 
is less obvious. Instead we approach the problem by systematically relating the Fock expan- 
sion to the matrix elements of a certain class of quark-gluon operators between the QCD 
vacuum |0) and the hadron states, taking advantage of the simplicity of |0) in light-cone 
quantization. The scale-dependence of the wave function amplitudes can then be traced to 
the wave function renormalization constants of quark and gluon fields. Following this, we 
derive the scale dependence of the parton densities from individual Fock components. The 
scale evolution of these contributions obeys an infinite set of coupled, linear differential- 
integral quations. When summed over all Fock contributions, we recover the well-known 
Dokshitzer-Gribov- Lipatov-Altarelli-Parisi (DGLAP) equation for parton densities. 

Before starting, let us remind the reader some salient features of light-cone quantization 
relevant for the following discussion [0]. The light-cone time x + and coordinate x~ are 
defined as x ± = l/y/2(x° ±x 3 ). Likewise we define Dirac matrices ^ = l/v^2(7° ±7 3 )- The 
projection operators for Dirac fields are defined as P± = (1/2)7 =F 7 ± . Any Dirac field ip can 
be decomposed into ip = ip + + tp_ with ip± = P±ip. ip + is a dynamical degrees of freedom 



and has the canonical expansion, 



+d[(k)v(k\)e i(k+ ^-^ ) ] . (1) 

Likewise, for the gluon fields in the light-cone gauge A + = 0, A± is dynamical and has the 
expansion, 

A 

+a[(k)e*(kX)e lik+ ^-^] . (2) 

■0_ and A~ are dependent variables. 

The key observation in this paper is that the light-cone Fock expansion of a hadron state 
is completely defined by the matrix elements of a special class equal light-cone time quark- 
gluon operators between the QCD vacuum and the hadron. These operators are specified 
as follows: Take the + component of the Dirac field ip+ an d the + _L-component of the 
gauge field F +± . [We sometimes label the _L components with index % = 1,2.] Assume all 
these fields are at light-cone time x + = 0, but otherwise with arbitrary dependence on other 
spacetime coordinates. Products of these fields with the right quantum numbers (spin, fla- 
vor, and color) define a set of operator basis. [This has been done in the past for light-cone 
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correlations in which all fields are separated along the light-cone Q.] Clearly, these opera- 
tors are not gauge-invariant because one cannot gauge-invariantize them by simply inserting 
string operators along the light-cone. Since all fields are at different points in the transverse 
directions, the operators do not have singularities requiring special renormalization. More- 
over, at equal light-cone time, there is no need to introduce a time-ordering among different 
fields because the difference is proportional to equal-light-cone-time commutators which are 
straightfoward to evaluate. In fact, to simplify the discussion, we assume all fields in the 
operators are normal ordered, i.e., the annihilation operators appearing at the right of the 
creations. We believe that the matrix elements of all these operators between the hadron 
state and the QCD vacuum yield complete information about the hadron wave function. 

As an example, let us consider 7r + meson with momentum P M along the ^-direction. The 
leading light-cone Fock states consist of a pair of up and anti-down quarks. The light-cone 
helicity of the ir + meson is zero, but the light-cone helicity of the quark-antiquark pair 
can either be zero or ±1. Use to represent the up-quark field in the coordinate 

space and d + (0) the anti-down-quark field. In the massless limit, the operator gL|_7 + 75M + 
yields a helicity-0 pair, and d + a +L ^^u + a helicity-1 pair. The helicity counting here is 
based on the chirality of the operators and the relation between chirality and helicity of 
massless fermions. [One can in principle construct similar operators without the 75 matrix; 
however parity forbids any finite matrix elements of them between the QCD vacuum and 
pseudo-scalar mesons.] The first operator defines a coordinate amplitude, 

(0|d + (0) 7 + 75« + (r,^)k + (P)) = o (r,£±)2P + , (3) 

where we normalize the state covariantly (P\P') = 2P + (2n) 3 5(P + - P + ')5 2 (P ± - P' ± ). 
Introducing the Fourier tranformation of the amplitude 

Mk±,x) = J ^re i(fc+r -^Vo(r,e±) , (4) 

we can invert Eq. (Q) to find 

W + (P)) = [ f 9 % - r^—M ^ fcO \bU x > £04i«(l - x, -k ± ) 

J {^) 2y/x{l — X) 1 

-bl^k^il-x,-^)] |0) + ... (5) 

where % is the color index. The creation and annihilation operators are normalized according 
to the commutation relation [b(k), b\k')} + = 2k + (2n) 3 5(k + - k + ')5 2 (k 1 _ - k' ± ). We have 
assumed here that the full QCD vacuum is a perturbative vacuum in light-cone quantization. 
In particular, we neglect the subtlety of zero-modes which might cause problems at some 
stage §. 

The operator with helicity-1 quark-anti-quark pair defines the amplitude 

(0|rf + (0) ( x +4 75« + (r,e±)|vr + (P)) = dVi(r,£±)2P + , (6) 

where i — 1, 2 is an index for transverse directions. Peforming a Fourier transformation on 
the both sides and inverting the equation, we find a light-cone Fock component 

{2n) 6 2a/x(1 - x) L 
+ (h + ifcO&y*, fc±)4^(l - ^ -k±)} |0> + ... (7) 
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The angular momentum content of the wave function is clear: For a quark-antiquark pair 
carrying helicity ±1, it couples to an orbital wave function with L z = =pl. Parity determines 
the relative sign of the two contributions. The phenomenological implications of (j)i(x,k±) 
for pion form factors have been discussed in the literature before ||. 

It is now straightfoward to study the cut-off dependence of 0o,i( x ; k±). For the moment, 
we focus on the transverse momentum cut-off A, although a cut-off in x is also needed at 
x — > in general. Besides the explicit cut-off dependence in the wave function amplitudes, 
the k± integration in Eq. (5) is implicitly bounded by A. In any cut-off scheme, the quark 
and gluon fields in QCD as well as the strong coupling constant a s depend on the cut-off. For 
large A, the dependence of quantum fields on A can be calculated in perturbation theory 
because of the asymptotic freedom. In fact, according to the standard renormalization 
theory, on has 

MO = 4 /2 (AM0, ^(0 = ZV\A)A"(£) (8) 

where and are indepedent of A, and Zf,a{A) are the wave function renormaliza- 

tion constants. Although the factorization in the above equation is scheme-dependent, but 
the A dependence itself is not. 

Going back to Eqs. (|3],|6]), it is now clear that the cut-off dependence of the wave function 
amplitudes <pi comes entirely from the field renormalization, 

<t>f(x,k ± ) =Z F (A)4i(x,k ± ) , (9) 

where <j>i(x, k±) is independent of A. 

We claim that the above feature holds for all components of the pion wave function in 
the Fock expansion. For instance, the most general two-quark, one-gluon wave function 
amplitudes can be defined through the matrix elements of the operators, 

d + (0)7 + 75i 71+J (^,r ?± ) M+ (r,e±) , 

d + (0)a +i l5 F +j (^, V± )u + (r,U) • (10) 
Their scale dependence comes entirely from the wave function renormalization constants 
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Z F (A)Z2 (A). In general, an n-particle Fock wave function amplitude with n q quark 
and antiquark and n g gluon creation operators has an explicit cut-off dependence through 
Z F q ^ 2 (A)Z^ 9 ^ 2 (A). Once again, all the momentum integrations are cut-off by A. 

Knowing the scale dependence of individual components of the hadron wave function, 
we can calcualte the scale dependence of their contributions to hadronic matrix elements. 
As an example, we consider in the remainder of the paper Feynman parton distributions, 
although the discussion applies to generalized parton distributions as well J7|. Deriving the 
parton evolution equation from light-cone wave functions has been considered in Ref. [BJ. 
Our approach here allows to uncover a set of new equations. 

Consider, for example, the two-particle wave- function contribution to the u quark distri- 
bution in the pion. We have 



u 2 (x,A) = 




Since at large k±, <fio{ x ,k±) g° es like l/k\ and <fii(x, kj_) like l/k\_ modulo logaritms [|TJ, 
the kj_ integration is convergent and hence A can be taken to infinity. Thus the only A 
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dependence in u 2 {x, A) comes from the wave function renormalization factor Zp. This yields 
the following evolution equation for U2 (x, A) 



;U2 (x,A) = -2 lF ^4^u 2 (x,A) , (12) 



dlnA 2 zv ' ' Ir 4tt 

where 7^ is the anomalous dimension of Zp and is gauge-dependent. [In physical gauges, it 
is positive-definite.] In light-cone gauge, 

lF = 2C F f dy 1 + {1 ~ y)2 , (13) 

Jo y 

where Cp = 4/3. The above integral diverges at y = ||, and we regulate the integral by 
cutting it off at y = e. Physics of course must be independent of any cut-off. 

Equation (12) indicates that the two-particle Fock state contribution to the up-quark 
distribution graduately diminishes as A — > 00. The physics is simple: As A gets larger, the 
state is probed is at shorter distance, and it becomes increasingly difficult for the meson to 
remain in the two-particle Fock component because of the radiation. As a consequence, the 
three-particle Fock amplitude increases at the leading-logarithmic level. In the light-cone 
gauge, "fp diverges at small x, and the radiation rate will depend on the cut off e. 

Consider now the three-particle Fock component contribution to the w-quark distribution. 
We write schematically, 

«3(*,A) = J ^^^^ \<f> A (x,k ± ,x',k' ± ,A)\ 2 , (14) 

where we have not considered individual quark-antiquark-gluon helicities and orbital angular 
momentum projections although this can be done straightforwardly. As discussed before, the 
wave function A (x, k±, x', k' ± ) has an explicit dependence on A through the wave function 

1 /2 

renormalization constant ZpZJ^ . Additional dependence comes from integrations over the 
transverse momenta k± and k' ± . 

We take derivative with respect to A using the chain rule. The derivative with respect to 
the wave function renormalization yields, 

; u 3 (x, A) = -(2 7F + lA )^lu 3 {x, A) + ... . (15) 



cilnA 2 ov ' ' y Ir 4vr 

where is the anomalous dimension of the gluon wave function renormalization. The 
physics of this part of the scale evolution is the same as the two-particle Fock component 
case: The splitting of the partons leads to the decrease of the probability for the pion to 
remain in the three-particle Fock state. 

The integrations over k± and k' ± do not yield divergences in general. In fact, there is no 
overall divergence (divergence arising from when all transverse momenta going to infinity 
at the same rate) because the power counting indicates that when two-momenta going 
to infinity at the same time, the integrals have negative superficial degree of divergence. 
However, there are sub divergences. These sub divergences arise from one-loop diagrams 
shown in Fig 1. The physics of these diagrams is that there are three-particle Fock amplitudes 
which are generated from the radiation of the two-particle Fock amplitude. Therefore, the 
result is proportional to m 2 (x, A). 
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FIG. 1: Contributions to the u-quark distribution from the three-particle component generated by 
the two-particle Fock component. 



First consider the loop intergral from the d-quark line shown in Fig. la. In this case, the 
k' ± integration is divergent whereas the k± integration is finite. It is easy to see that when 
A changes, the result 3-particle distributions also changes. From this diagram, one finds, 



d 



din A 2 



, 4X oc s (A) , . . 

u 3 {x, A) = ... + If—. u 2 {x, A) . 

An 



(16) 



where the plus sign indicates that the three-particle Fock component receives a contribution 
from the radiation of the two-parton states. 

Finally, let us consider the gluon radition from for the u-quark line as shown in Fig. 
lb. The integration over kj_ is now divergent whereas the one over k' ± is finite. The k± 
integration can be done using the standard light-cone perturbation theory, 



d 



din A 2 



u 3 (x, A) 



a s (A) [ 1 dy l + y 2 
■■■ H ^-C F / u 2 



2tt 



y i -y 



X A 

-A 

y 



(17) 



where the divergence at y — 1 must be regulated. Adding everything together, we obtain 
the complete evolution equation for u 3 is 



d 



din A 2 



u 3 (x,A) 



q B (A) 

Air 
+2C F 



- (2'j F + 7a)«3(^, A) + 'j F u 2 (x, A) 
1 dy 1 + y 2 fx 



y 1-2/ 



u 2 



A 



y 



(18) 



This is an inhomogeneous equation with a driving term u 2 . 

Going to wave function amplitudes with four and more partons posts no special difficulty, 
except one has to take into account the mixing with the singlet contribution. Take the 
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example of four-parton amplitudes for which three flavor structures udgg, uduu, and uddd 
must be considered separately. For udgg, the gluon radition from u, d an g of the three- 
parton component udg yields, 



d 



din A 2 



u 



udgg 



(x,A) = 



<* a (A) 

47T 



-2( 7 f + 7a)< m (x,A) 



+ (7f + 7aiK(*, A) + 2C F f d JL\±f 

Jx y i- 



-« 3 



A 



2/ 



(19) 



where 7^1 is the part of the gluon anomalous dimension from the gluon loop. On the other 
hand, for uduu, the gluon splitting into uu pair yields, 



Q= 8 (A) 

47T 



-4 7F <" u (x,A)+7 A2M3 (x, A) 



+2 f d ^T F {y^ + {I - yf)g, (*,A 

Jx y \y 



(20) 

where Tp = 1/2 and 7,42 is the part of the gluon anomalous dimension from the uu loop. 
g 3 (x,A) is the gluon distribution from the udg Fock amplitude. Finally, for uddd, the gluon 
splitting into dd pair yields, 



d 



dlnA 2 " 4 



uddd. 



x,A) = 



a s (A) 
An 



A lF uf dd {x,A)+ lA3 u 3 {x,A) 



(21) 



where 7^3 is the part of the gluon anomalous dimension from the dd loop. When adding all 
the contributions (7^ = 7ai + lA2 + Iaz) , we have, 

a s {A) 



d 

din A 2 



u±(x, A) 



47T 

+2 



- 4 7F <^(x, A) - 2( 7f + 1a)uT"(x, A) + ( 7F + lA )u 3 {x, A) 



1/ I 



L-A- «3 

1-2/ 



y 



,A)+T F (y' + (l-yY)g 3 (^,A 



x 



(22) 



To keep the evolution simple, we also consider the anti-up quark distribution at this 
order. The only contribution is from uduu component for which we have 



d _ , . N oc s (A) 
dh^ U4{x > A) = — 



-47 F U 4 (x, A) + 2 



dy 

y 



T F (y' + (l- y y)g 3 (-,A 



x 



y 



(23) 



Therefore, if we define the valence up-quark distribution u v = u — u, then 



d , . a s (A) 
u Av {x,A) = — 

47T 



din A 2 



-A lF ut q ~ q {x, A) - 2{ lF + lA )ut" (x, A) 



udgg. 



+ {1f + 1a)u 3 {x,A) + 2 [ ^l±^U 3 f-,A 

Jx y i-y V?/ 



(24) 



without the complication from the g — > qq kernel. 

It is not difficult to see that the evolution equation for u nv from Fock states with n partons 

is 

n n—2 



d 



din A 2 



U n v(x, A) 



a s (A) 
2C F 



- ^2 7iUnv(x, A) +^ liU n ~lv(x, A) + 
i=l 

dyl + y 2 



i=i 



, x 



-u n - lv (-,A) 
y 1-2/ 2/ 



(25) 
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where the sum over n — 2 ji excludes one 7^ and one 7^. The first term in Eq. (25) describes 
the depletion of the n-particle Fock component due to the gluon radiation into n + 1 -particle 
component; the second and third terms describe the increase of the n-particle component 
due to the gluon emission of the n — 1-particle component. The difference between the later 
two comes from whether the gluon is radiated from the active particle or the spectators. 
The total u v (x) distribution is a sum over all possible Fock components, 

00 

u v(x) = y\ fe Qr) • (26) 

8=2 

Summing over all the equations for the individual Fock components, we recover the standard 
DGLAP equation, 

d i A x a 8 (A) f 1 dy 1+y 2 fx \ 

; u v {x,A) = — — / — - — u v -,A (27) 



dlnA 2 " 2tt J x y (1 - y) + \y 

which is an important check. Using the same procedure, one can derive evolution equations 
for other types of parton distributions, such as the singlet quark and gluon distributions, 
quark helicity and transversity distributions, as well as higher-twist distributions. 

In summary, we find that the light-cone wave functions of hadrons in QCD can be entirely 
determined by the matrix elements of a class of quark-gluon operators. From this, we derive 
an infinite set of evolution equations for the parton distributions contributed by individual 
n-particle Fock components. These equations are consistent with the well-known DGLAP 
equation and are useful for phenomenological studies of hadronic structures and model 
buildings. They should also provide important constraints for wave functions derived from 
light-cone quantization. 

Note added: After this paper was finished, we were informed that Ref. [jTO] has studied the 
evolution of the 2 and 3 particle Fock components, using an explicit light-front calculation. 
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